Tomography of the atmosphere of long-period variable stars 
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Abstract. This paper reviews the results obtained so far with the newly developed tomographic technique, 
which probes the velocity field across the atmosphere of long-period variable (LPV) stars. The method 
cross-correlates the optical spectrum with numerical masks constructed from synthetic spectra and probing 
layers of increasing depths. This technique reveals that the line doubling often observed in LPV stars around 
maximum light is the signature of the shock wave propagating in the atmosphere of these pulsating stars. 

1 Introduction 

It is known since long that the brightness variations of long-period variable stars (LPVs) go along 
with spectral changes such as the doubling of the absorption lines around maximum light [1,3,10]. 
This line-doubling phenomenon appears to be shared by all LPVs [5] , and occurs when the shock 
wave associated with the envelope pulsation propagates through the line- forming region [3,17]. 
Since it is almost impossible to study individual line profiles in the very crowded optical spectra of 
LPVs, it is necessary to resort to a technique that correlates the observed spectrum with a spectral 
template [4,6,16]. The shape of the resulting cross-correlation function (CCF) thus represents 
the average shape of those lines in the observed spectrum that match the template. By carefully 
designing spectral templates collecting all spectral lines forming at given depths in the atmosphere, 
it is even possible to probe the velocity field in pulsating variables. This is the general principle 
guiding the tomographic 1 method described in this review. When applied to a temporal sequence 
of spectra, the method is even able to reveal in a spectacular way the outward propagation of 
the shock wave, as shown in Sect. 4 and Figs. 5-6! 

2 Tomography of the atmosphere 

The study of the velocity field in the atmosphere of LPVs poses a special challenge, as the spec- 
trum of these stars is extremely crowded, particularly in the optical domain. The cross-correlation 
technique provides a powerful tool to overcome this difficulty. The information relating to the line 
doubling (velocity shift and line shape) is in fact distributed among a large number of spectral 
lines, and can be summed up into an average profile, or more precisely into a cross-correlation 
function (CCF). If the correlation of the stellar spectrum with a mask involves many lines, it is 
possible to extract the relevant information from very crowded and/or low signal-to-noise spectra. 
The CCF writes 

CCF(AX) = [ 2 s(A - AX) m(A) dA, (1) 

where s(A) is the observed spectrum, m(A) is the template (a binary template has been adopted, 
being around spectral lines of interest, and 1 elsewhere), and Ai and A2 are the boundaries of 

1 The word tomography is used here in its etymological sense ( 'display cuts ') , which differs somewhat 
from the broader sense in use within the astronomical community (reconstruction of a structure using 
projections taken under different angles). 
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the spectral range covered by the observed spectrum. The radial velocity V r is then obtained from 
the wavelength shift Z\A(min CCF) where the CCF is minimum: 

ZU(min CCF) 

Vr = 0-5 (A 1+ A 2 ) C ' (2) 

where c is the speed of light in vacuum. We refer the reader to reference [6] for a detailed description 
of the CCF mathematical properties. 

The tomographic method rests on our ability to construct reliable synthetic spectra of late-type 
giant stars [13-15], and from those, to identify the depth of formation of any given spectral line. 
Rigorously, the contribution function (CF) to the flux depression [2] should be used to evaluate 
the geometrical depth at which a line forms. However, it would be a formidable task to compute 
the CF for each line appearing in the optical spectrum of LPV stars. For the sake of simplicity, the 
'depth function' x = x(X) is used instead, which provides the geometrical depth corresponding to 
monochromatic optical depth t\ = 2/3 at the considered wavelength A. This function expresses the 
depth from which the emergent flux arises, in the Eddington-Barbier approximation, and should 
not differ much from the average depth of formation for sufficiently strong lines [11]. 

Different masks Ci are then constructed from the collection of Ni lines \ij(l < j < N) such 
that Xi < x(Xij) < CGj+i = Xi + Ax, where Ax is some constant optimized to keep enough lines 
in any given mask without losing too much resolution in terms of geometrical depth. A natural 
limit to the resolution that can be achieved on the geometrical depth is provided by the width of 
the CF [4] , and according to the Nyquist-Shannon theorem of elementary signal theory, there is no 
advantage in taking Ax smaller than half the CF width of typical lines probed by the mask. Each 
mask d should then probe lines forming at (geometrical) depths in the range Xi , Xi + Ax in the 
atmosphere. These masks are used as templates to correlate with the observed spectra of the LPV 
stars. This procedure should provide the velocity field as a function of depth in the atmosphere 
(Fig. 1), whereas for static atmospheres, all CCFs should yield the same radial velocity (Fig. 3). 

A more detailed description of the method can be found in [4] , while specific tomographic masks 
are available at http://www-astro.ulb.ac.be/Html/home. html#tomography. The set of tomographic 
masks used in this paper was constructed from a synthetic spectrum at T c g — 3500 K and log g = 0.9 
(see reference [4] for details). The properties of the masks (in terms of numbers of lines and depths 
probed) are listed in Table 1. The depth x is in fact expressed in terms of a reference optical 
depth x = log t , where To is the optical depth at the reference wavelength of 1.2 /im. There is a 
one-to-one correspondence between log To and the geometrical depth. 



Table 1. The synthetic templates constructed from the model at 3500 K and logg = 0.9. 



Mask # 


depth probed 


# of lines per mask 




X = log T1.2 M m 




c\ 


Xi,Xi + Ax 


Ni 


Ci (innermost) 


-2.00,-2.75 


777 


c 2 


-2.75,-3.50 


610 


c 3 


-3.50, -4.25 


433 


c 4 


-4.25, -5.00 


321 


c 5 


-5.00, -5.75 


168 


c 6 


-5.75, -6.50 


167 


cv 


-6.50, -7.25 


94 


Cs (outermost) 


-7.25, -8.00 


46 
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Fig. 1. The CCFs of the Mira Z Oph at phase 0.08 obtained with the tomographic masks. Note how the 
shape of the CCFs evolve from the innermost layer (involving ascending matter only, hence CI exhibits a 
single blue peak) to the outermost layer (involving mostly matter falling in, hence C8 exibits predominantly 
a red peak). This spatial sequence of line doubling reflects the presence of a shock wave in the line- forming 
region, with the shock front being centered on the layer probed by the mask C5 




Fig. 2. Comparison of the A 6358.69 Fel line as seen in the spectrum of Z Oph at phases 0.08, 0.17, 1.04 
and 1.21 (bottom row, from left to right), with the CCF obtained with the tomographic mask C5 (upper 
row), to which belongs the A 6358.69 Fel line. Z Oph is warm enough at maximum light (K3) for its 
spectrum to be not too crowded so that clean, almost unblended, lines as the one displayed here may be 
isolated 



An a posteriori validation of the method is provided by the fact that (i) the CCFs obtained 
with the tomographic masks for the non-pulsating M giant fi Gem (Lb variable with a variability 
range of only 0.11 mag in the Hipparcos catalogue) are all single-lined and yield almost identical 
radial velocities, as expected (Fig. 3), and (ii) the asymmetric CCFs obtained for the warm Mira 
variable Z Oph reflect the shape of the individual lines probed by any given mask (Fig. 2). The 
asymmetric, and sometime even double-peaked, CCFs observed for Mira variables thus correctly 
represent the underlying spectral lines, and are not artifacts of the correlation process. Moreover, 
around maximum light, the shape of CCFs in the successive tomographic masks follows a definite 
sequence: the CCF in the deepest layer exhibits a single blue component whereas a single red 
component is observed in the outermost CCF, and double-peaked CCFs are observed in between 
(Fig. 1). This sequence, which is observed in almost all Mira around maximum light, is consistent 
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Fig. 3. Same as Fig. 1 for the non-Mira star /j, Gem. Note how the CCFs broaden and become more 
contrasted in the upper atmospheric layers (C8 is the outermost mask), reflecting the fact that the strongest 
spectral lines form in the outermost layers. The radial velocity (in km/s) derived from each CCF is given 
in each panel. 




Fig. 4. The Schwarzschild scenario: temporal sequence followed by the intensity of the red and blue compo- 
nents of absorption lines close to maximum light, when the shock wave propagates through the photosphere 



with the Schwarzschild scenario describing the effects of the passage of a shock wave through the 
atmosphere of a pulsating star (Sect. 3). 

3 The Schwarzschild scenario 

The passage of a shock wave through a stellar atmosphere is expected to give rise to a very 
specific temporal evolution of the spectral line shape (Fig. 4), as first noted by Schwarzschild 
[17] in relation with W Vir Cepheids. When the shock wave penetrates the line-forming region, 
the velocity discontinuity associated with the shock gives rise to the doubling of the spectral 
line: the red component of the line is formed in the matter falling in above the shock, whereas 
the blue component arises from the ascending matter lying behind the shock. The intensity of 
the blue and red components of a double line observed around maximum light in a LPV star 
should follow the temporal sequence illustrated in Fig. 4. This is indeed the case, as revealed by 
the evolution of the CCFs of the Mira variables RT Cyg and RY Cep around maximum light 
(Sect. 4). But the tomographic method described in Sect. 2 allows us to go one step further, by 
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showing that the Schwarzschild scenario also reveals itself in terms of a definite spatial sequence of 
profiles (Fig. 1). By combining those spatial sequences obtained at successive phases in the light 
cycle, the tomographic method reveals that line-doubling appears later in upper layers of the Mira 
atmosphere, thus translating the upward motion of the shock wave (Sect. 4). 

4 Application to the Mira variables RT Cyg and RY Cep 

A long-term monitoring of the Mira stars RTCyg (P = 190 d; 6.0 < V < 13.1) and RY Cep 
(P = 145 d; 8.6 < V < 13.6) has been performed with the fibre-fed echelle spectrograph ELODIE 
[7]. The spectrograph ELODIE is mounted on the 1.93-m telescope of the Observatoire de Haute 
Provence (France), and covers the full range from 3906 A to 6811 A in one exposure at a resolving 
power of 42 000. 

For RT Cyg, a monitoring in August-September 1999 covered phases —0.20 to 0.16 around 
maximum light, with 32 spectra obtained during this phase range (corresponding to an average 
resolution of A(f> = 0.01; Fig. 5). For RY Cep, the monitoring covered a full light cycle extending 
from August 2001 (phase -0.17) to February 2002 (phase 0.92), and 40 spectra were obtained 
(Fig. 6). 

The sequences of CCFs presented in Figs. 5 and 6 clearly obey the Schwarzschild temporal 
evolution, with a single red component transforming progressively into a single blue component 
around maximum light. There are moreover clear phase lags between the different layers, this 
transformation occurring at later phases in outer layers. Thus, the Schwarzschild scenario holds for 
Mira variables, and this conclusion definitely points towards the velocity stratification associated 
with the shock wave as the cause of the double absorption lines observed in Mira variables, as 
opposed to complex radiative processes (e.g., radiative release of thermal energy into the post- 
shock layer or temperature inversion) [8,9]. 

Radial velocities have been extracted from the CCFs of Figs. 5 and 6 (top panels) to yield the 
curves displayed in the bottom panels of the same figures. We defer the detailed comparison between 
the radial- velocity curves of RY Cep and RT Cyg, and their interpretation, to a forthcoming paper. 
We just stress here noteworthy features on the radial-velocity curve of RY Cep: (i) the red peaks 
(corresponding to matter falling in) are observed at the same velocity in all three masks; they 
disappear at increasingly later phases as one considers layers closer and closer to the surface, since 
the outermost layers are the last ones to be penetrated by the shock wave which suppresses the 
infalling component (i.e., red peak); (ii) the maximum outward velocity is the same for the three 
masks (~ —247 km/s), and is reached around phase —0.1. After that, matter will decelerate in the 
innermost layers first, as they lose the impetus provided by the shock wave which is moving away. 

5 Conclusion and perspectives 

The tomographic technique presented in this paper opens new perspectives for the study of the 
dynamics of LPV stars, as it allows to visualize the outward motion of the shock wave in the 
atmosphere. At this stage, however, the results remain qualitative as it is difficult to quantify the 
geometrical depth of the layers probed by the different masks. It is hoped that, in the future, 
this method may be combined with interferometric imaging in the different masks to derive the 
corresponding radii. This would give direct access to the velocity of the shock, for instance. 
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Fig. 5. Top panel: Sequence of cross-correlation profiles of RTCyg obtained with the tomographic masks 
C3, C4, C5 and C7 (from left to right) in August-September 1999, around maximum light. The labels 
beside each CCF refer to the corresponding visual phase based on the AAVSO ephemeris [12] 
Bottom panels: Radial velocity curves obtained with the tomographic masks C3 (leftmost panel - solid line 

- triangles), C4 (middle panel - long-dashed line - squares), and C7 (rightmost panel - short-dashed line 

- circles). The C4 radial- velocity curve (long-dashed line) has been duplicated in the left- and rightmost 
panels to allow an easy intercomparison 




Fig. 6. Top panel: Same as Fig. 5 for RY Cep with the tomographic masks C3 (inner), C4, C5 and C6 
(outer), from left to right. The labels beside each CCF refer to the corresponding visual phase based on 
the AFOEV monitoring. 

Bottom panels: Radial velocity curves obtained with the tomographic masks C3 (leftmost panel - solid line 

- triangles), C4 (middle panel - long-dashed line - squares), and C6 (rightmost panel - short-dashed line 

- circles). The C4 radial- velocity curve (long-dashed line) has been duplicated in the left- and rightmost 
panels to allow an easy intercomparison 
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